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Abstract 

We present 3D models of the Galactic magnetic field including regular 
and turbulent components, and of the distribution of matter in the Galaxy 
including relativistic electrons and dust grains. By integrating along the 
line of sight, we construct maps of the polarized Galactic synchrotron and 
thermal dust emissions for each of these models. We perform a likelihood 
analysis to compare the maps of the Ka, Q, V and W bands of the Wilkinson 
Microwave Anisotropy Probe (Wmap) and the 353 GHz Archeops data to 
the models obtained by varying the pitch angle of the regular magnetic field, 
the relative amplitude of the turbulent magnetic field and the extrapolation 
spectral indices of the synchrotron and thermal dust emissions. The best-fit 
parameters obtained for the different frequency bands are very similar and 
globally the data seem to favor a negligible isotropic turbulent magnetic field 
component at large angular scales (an anisotropic line-of-sight ordered com- 
ponent can not be studied using these data). From this study, we conclude 
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that we are able to propose a consistent model of the polarized diffuse Galac- 
tic synchrotron and thermal dust emissions in the frequency range from 33 
to 353 GHz, where most of the CMB studies are performed and where we 
expect a mixture of these two main foreground emissions. This model can be 
very helpful to estimate the contamination by foregrounds of the polarized 
CMB anisotropics, for experiments like the Planck satellite. 
Keywords: ISM: general - Methods: data analysis - Cosmology: 
observations - millimeter 



1. Introduction 



The Planc Kill satellite mission ( iPlanck- Collab or at ioni ( l2011al ) 



Tauber et al, 



(2010|)), currently in flight, will provide measurements of the CMB anisotropics 



both in temperat ure and polarization over t he full sky with an unprece 



20051). In particular, it should be 



dented accuracy (iPlanck-CoUaborationl . 
able to m easur e the t e nsor-to-scala r ratio , r, which sets the energy scale of 
inflation Jhvthl (|l98J); 



Peiris et al. 



(120031 )). For a extended mission and in 



the c ase of no direct dete c tion, Planck should se t an upper limit of r < 



0.03 ( iBurigana et al. 



2010 



Efstathiou et al. 



2009 



Efstathiou and Grattonl. 



Baumann 



20091)) . and thus provides tight constraints on inflationary models ( 
( I2OO9I )). To achieve this high level of sensitivity, it is necessary to estimate 
accurately the temperature and polarization anisotropics from foreground dif- 
fuse Galactic emissions and from point-like and compact sources of Galactic 
and extraCalactic origin. A reliable estimation of the residual contamination 
due to the foreground emissions after application of component separation 



"'^http://www.rssd.esa.iiit/index.php?project=Planck 
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methods (see iBetoule et al.l ( l2009l ): iLeach et al.l (l2008[ l for recent studies) is 
thus necessary to retrieve the cosmological information from the Planck 
data. 



As summarized in lFauvet et al.l ( 120111 ) the main polarized fore ground con 



tribu tions come from the diffuse Galacti c synchrotron emission (IPage et al. 



2007 ) and from thermal dust emission (IBenoit et al. 



2004 



Ponthieu et al. 



20051). The polarize d synchrotr o n ern is sion have a l ready been modelled 



Page et al 



(120091 



2004 



20071) : 



Fauvet et al 



Sun et al 



(|2008[ ): 



Jaffe et al 



mm; 



Jansson et al 



(120111 ) based on models of the Galactic magnetic field ( iHan et al. 



20061) and of the relativistic electrons in the Galaxy. Concerning the 



thermal dust emission, modeling based on the physical or i gins o f this polar- 



i zed emission has already be en discussed in 



Ponthieu et al. 



(I2007h 



Fauvet et al. 



(120 111 ). 



(120051 ) : IPage et al. 



the 3D ioint r node 



Fauvet et al 



of the Galac- 



(1201 If ). In the pre- 



We propose here an extended analysis o 
tic polarized diffuse emissions discussed in 
vious analysis we focused on the Wmap data at 23 GHz and Archeops 
at 353 GHz data where synchrotron and thermal dust emissions dominate, 
respectively. Here, we use complementary data : the other frequency bands 
of Wmap from 33 to 94 GHz, where a mixing of those emissions is expected, 
and the Archeops data at 353 GHz. Furthermore we apply here a pixel-to- 
pixel likelihood based compa r ison i nstead of a Galactic profile-based method 



as discussed in 



Fauvet et al. 



( 1201 ih 
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The paper is structured as follows: Sect. [2] describes the five-year Wmap 
and Archeops data set used in the analysis. In Sect J3] we describe in detail 
models for the polarized foreground emissions. Section H] discusses the 3D 
model of the Galaxy using to construct the polarized Galactic emissions. The 
models are statistically compared to the data in Section Eland we discuss the 
results in Section [61 We finally conclude in Section [71 

2. Observations 

2.1. Diffuse Galactic synchrotron emission 

The synchrotron emission is an important contributor to the diffuse sky 
emission at both radio and microwave observation frequencies. 



In intensity, the 408MHz all-sky continuum survey (IHaslam et al.l ( 1l982l )). 
at a resolution of of 0.85 degrees, is a good tracer of the synchrotron emission 
and it will be used in the following as a template. In particular, we use the 



408MHz all-sky ma p available on the 



pixelisation scheme (iGorski et al 



;AMBDA website in the HEALpix 



20051 ). We correct this map from the con- 
tribution from free-free emission which estimated to be about 30 %. The 
corrected maps was then downgraded to the resolution of the WMAP and 
Archeops maps discussed below. 



In polarization, Faraday rotation introduces complications into the inter- 
pretation of the radio data since strong depolarization is observed for frequen- 



cies lower tha n 10 GHz, p a rticu 



flBurnl mm -. 



Sun et al 



mm 



ar 



y concerning the in ner part of the Galax y 



Jaffe et al. 



mm-. 



Jansson et al 



(120091); 



La Porta et al.l ( 120061 )). A detailed discussion on the depolarization process, 



basec. 



on data from the Leiden survey, can be found in 



La Porta and Burigana 



(120061 ) ■ For this reason the best polarized Galactic d iffuse synch r otron trac- 



ers are at high frequency such as the Wmap survey (jPage et al.l (120071 )) 



We used here the five-year Wmap Q and U low resolution maps for the 
frequency channels centered at 3 3 (Ka band), 41 ('Q band ), 61 (V band) and 



Gold et al. 



(120091 ) ). These data are 



94 GHz (W band) (iPage et al.l (120031); 
available on the LAMBDA website in the HEALpix pixelisation scheme at 
Nside = 16. The associated noise is estimated using the full noise correlation 
matrix also available on the LAMBDA website in the same resolution. 

2.2. Thermal dust 

The thermal dust emission is significant in the Wmap data only for fre- 
quencies above 70 GHz, then we also used here the Archeops 353 GHz Q 
and U maps as tracers of th e polarized thermal dust emission. Those maps 
cover about 20 % of the sky (IMacias-Perez et al.l . 120071 ) and were filtered and 
downgraded to Nside = 16 to make them comparable to the Wmap ones. 



In intensity, the most accurate me asurements of the thermal dust emis- 
sion are those of the IRAS satellite ( iNeugebauer et al.l . Il984j ) and in par- 
ticular at 100 fim. We use here predicted full-sky maps of sub millime- 
ter a nd microwave em i ssion from the diffuse interstellar dust in the Galaxy 



from 



Finkbeiner et al. 



( 119991 ) which were produced cor abining the IRAS dat a 



at 6.1 arcmin and the COBE DIRBE data at 40 arcmin (ISchlegel et al. 



1998h 



These maps were downgraded to the resolution of the WMAP and Archeops 



maps presented above. 



3. Emissivity model in polarization 

We present in this section a realistic model of the diffuse polarized syn- 
chrotron and dust emissions using a 3D model of the Galactic magnetic field 
and of the matter density in the Galaxy. We will consider the distribution 
of relativistic cosmic-ray electrons (CREs), ucre, for the synchrotron emis- 



sion and t 
Following 



le distribution of du st grains, ridust, for the thermal dust emission. 

( 120111 ) we calculate the Stokes parameters I, Q and 



Fauvet et al 



U for the Galactic polarized emission along the line of sight as follows. 



For the synchrotron emission (IRybicki and LightmanI (119791 )) we write 

1/ 



/r^(n) 



n 



0.408 

V 

0.408 



^Has/ff(n) 

/ cos(27(n, s))ps (-BF(n, s) + 5|(n, s)) ncRE(n, s)ds 



c/^y°^(n) 



/i 



Has/ff 



/ (Sf (n, z) + B2(n, s)) ncRE(n, s)ds 



(n) 



0.408 



/ sin(27(n, s))p, (-B2(n, s) + 5f (n, s)) ncRE(n, s)ds 



(1) 
(2) 

(3) 

(4) 

(5) 



/ (Bf (n, s) + S2(n, s)) ncRE(n, s)ds 
where i?„(n, s) is the magnetic component along the line-of-sight n, and 

-B;(n, s) and i?t(n, s) the magnetic field components on a plane perpen- 
dicular to the line-of-sight. Notice that the 3 vectors n,l,t form an or- 
thonormal basis being / and i oriented to the north and to the east respec- 
tively in_a^_glane_2erpendku n. The polarization fraction is set to 
75% (IRybicki and LightmanI . Il979l ). The polarization angle 7(n, s) is given 
by : 
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7(n, s) 



- arctan 



2Biin, s)-Btin, s) 



(6) 



The distribution of relativistic electrons, ncRE, is described in detail in 
section |H luas/s is the reference r nap in intensi t y con structed from the 



4O8 MHz all sky continuum survey 



Haslam et al. 



fll982[ ) after subtraction 



of the bremsstrahlung (free-free) emission and u is the frequency of ob- 
servation. To subtract the free-free contribution we used the Wmap K- 



band fre e-free foregroun c 



(MEM) flHinshaw et al. 



map 



2007 



generated from the m aximum entropy method 



Bennett et al. 



2OO3I ). Notice that we do no t 



use the synchrotron MEM intensity map at 23 GHz (IHinshaw et al.l (120071 )) 
as a synchrotron template to avoid any possible Anomalous Microwave Emis- 
sion (AME) contamination (the Wmap team made no attempt to fit for the 
latter). The spectral index /J^ used to extrapolate maps at various frequencies 
is a free parameter of the model. The SED of the synchrotron emission 
in the radio and microwave domain, and in particular in the 33 to 
353 GHz range, can be w ell approximated by a p ower law in an- 



tenna temperature units (iRybicki and LightmanI (119791 )). This is 



due to the fact that the energy spectrum of the Galactic relativistic 
electrons producing the radio and microw ave synchrotron e missio n 
is also well approximated by a power law (jKobayashi et al.l (l2004j )). 
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For the thermal dust emission we write 



Q 



dust / 



n 



Ifds 



Ifds 



' ' V353/ 
' ' V353/ 



/ cos(27(n, s)) sin^(a)/„ormPdndust(n, s)ds 



/ ndust(n, s)ds 



u: 



dust / 



n 



4ds 



V3537 



/ sin(27(n, s)) sin^(Q;)/normPdndust(n, 
/ ndust(n, 



(7) 
(8) 
(9) 
(10) 

(11) 



where the dust polarization fraction pd is set to 10 % (jPonthieu et all l2005l ) 
based on the Archeops data, and Udustl'";^) is the dust grain distribution 
discussed in section |H The sin2{a) term accounts for the geometrical sup- 
pression and fnorm is an empirical factor which accounts for the misalign- 



ment between dust grains and the magnetic field lines (see (iFauvet et al 



20111) for details ) . Th e reference map, Ifds was taken to be model 8 in 



Finkbeiner et al. 



(119991 ) at 545 GHz. The spectral index 13d used to extrap- 
olate maps at various frequencies is a free parameter of the model. In the 
following we work on antenna temperature, Rayleigh- Jeans units. Assuming 
nearly constant dust temperature across the sky then a power-law approxi- 
matio n for the thermal dus t emiss ion in antenna temperature units can be 



used (jPlanck- Collab or at ion 



2011c). 



4. A 3D modeling of the Galaxy 

We describe here the 3D model of the Galaxy as used in the previous 
Stokes parameter definitions both for synchrotron and dust. 
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4.1. Matter density model 

In galactocentric cylindrical coordinates (r, z, (p) we consider an exponen- 



tia 
by 



distribution of relativ i stic e lectrons ucre on the Galactic disk motivated 



Drimmel and Spergell (120011 ): 



ncRE{r,z) 



(12) 



cosh2{z/ncRE,h)' 

where ncRE,h defines the width of the distribution vertically and it is set to 1 
kpc in the following. ncRF.r - defines the distrib ut ion radially and it i s set to 



3 kpc (see ISun et all ( 12008 



Jaffe et al. 



value of no,e is set to 6.4 x 10 ^.cm ^ ( ISun et al 



torn 



Fauvet et al 



20081 ). 



(l201l[ )l The 



The density distribution of dust grains in the Galaxy is poorly known and we 
therefore choose to describe it in the same way as for relativistic electrons: 



nrf(r, z) = no. 



d ■ 



(13) 



cosh.2{z /nd,h) 

where nd,r and nd,h are the radial and vertical widths of the distribution. 
In the following we set them to 3 and 1 kpc respectively. Notice that we 
expect these two parameters to be strongly correlated for both dust grain 
and electr on distribu t ions, h ence we deci d ed to fix one of them as in previous 



analyses (ISun et al. 



(120081) 



Jaffe et al 



(120101 )1 We have tested different 



values of these two parameters and found no impact on the final results. 
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4-2. Galactic magnetic field model 

The Galactic magnetic field model consists of a regular component and a 
turbulent component such that 



B 



toil 



Breg(l") + ^turb ^turbi''^) 



(14) 



where Afy^j-b is the amplitude of the turbulent component and it is a 
free parameter of the model. In the following we will express it as 
a relative amplitude with respect to the amplitude of the regular 
component. 



4-2.1. Regular component 

The regular part Breg(r) is a M odified Logarithmic Spiral model, discussed 



in detail in 



Fauvet et al. 



teOllI ). In galactocentric cylindrical coordinates 



[r, 0, z) it reads 



B(r) = Bregir)[cos{(j) + (3) In ^-^j sin(p) cos(x) • 

— cos(0 + /3) In I — ) cos(p) cos(x) ■ 
+sin{x) ■ Uz], 



(15) 



where p is the pitch angle, /3 = 1/ tan(p) and tq is the radial scale set to 
7.1 kpc. x(r) = \r,(r)(z/ zr,) is the vertical s c ale, w ith xo = 22.4 degrees 



and zq = 1 kpc. Following 



Taylor and Cordes 



( I1993I ) we restrict our model 



to the range 3 < r < 20 kpc. The lower limit is set to avoid the center 
of the Galaxy for which the physics is poorly constrained and the model 
diverges. The intensity of the regular field is fixed using pulsar Faraday 
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rotation measurements by 



HanetaL 



(120061 ) 



Breg{r) 



Bo e «s 



(16) 



where the large-scale field intensity at the Sun position is -Bq = 2.1 ± 



0.3/iG and the physical cut Rb = 8.5 ± 4.7 kpc. 



Sun and the Galacti c cente r, Rq is set to 8 kpc (lEisenhauer et al. 



'he distance between the 



fl2003h : 



Reid and Brunthalei] (120051 )1 



4-2.2. Turbulent component 

In addition to the large-scale Galactic magnetic field, Faraday rotation 
measurements on pulsars in our vicinity have revealed a tu r bulen t component 
on scales smaller than a few hundred pc flLvne and SmithI (119891) '). Moreover 



it seems to be present on large angular scales ( iHan et al 



(120041 )) with an am- 



plitu de estimated to b e of the same order of magnitude as that of the regular 
one (IHan et all (|2006[ )). The magnetic energy in Fourier space, Esik), asso- 
ciated with the turbulent componen t is well described by a power spectrum 



of the form (IHan et all (12004 



20061)) 
Esik) = C 



(17) 



where a = —0.37 and C = (6.8 ±0.3) • 10~^'^ergcm~^kpc. To obtain the 3D 
spatial distribution of the turbulent magnetic field we produced independent 
Gaussian simulations from the above power spectrum in the x, y and z di- 
rections on boxes of of 512^ pixels at a resolution of 56 pc. We renormalize 
the simulated boxes so that the total amplitude of the turbulent component 

is Aturb- 
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Notice that we do not includ e here an anisotropic /orde red component as 



discussed in 



Jaffe et al 



torn . As in flFauvet et al 



2011 



), o ur reg ular com- 



Jaffe et al 



torn call the 



ponent is then equivalent to the sum of what 
coherent and ordered fields. The latter, also called ordered random compo- 
nent, will be not considered in this paper because it can not be distinguished 



from the large-scale magnetic field when studying polarization intensity only. 



5. Method 



We compute I, Q and U maps for the synchrotron an d thermal dust 



emissi ons with a modified version of the Hammurabi code (jWaelkens et al 



(120091 ) ). Each map is generated 



jy integrating in 10 steps along each line-of- 



sight defined by the HEALPlx ljGorski et al. 



20051 ) iVgide = 16 pixel centres. 



The integration continues out to 25 kpc from the observer situated 8.5 kpc 
from the Galactic centre. These full-sky maps are computed for a grid of 
models obtained by varying the pitch angle, p, the turbulent component am- 
plitude, Aturb and the spectral indices of the synchrotron and thermal dust 
emissions Ps and (3d. The latter are assumed to be spatially co nstant on the 
sky. Dealing with a inore re alistic varying spectral index (see 



fl2007h 



La Porta et al. 



Kogut et al 



(120081 ) for detailed studies) is beyond the scope of this 
paper. However, we ensured that this hypothesis does not impact the results 
for the other free parameters in the model. Indeed, we produced simulated 
Wmap observations with spatially varying synchrotron spectral index and 
analyzed them assuming a constant one. No significant bias was observed 



for any of the other parameters and the error bars were equivalent to those 
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Table 1: Parameters of the 3D Galactic diffuse emissions model. 



Parameter 


Range 


Binning 


P (deg) 


[-80.0,15.0] 


5.0 




[-4.5,-2.8] 


0.05 


(3d 


[0.05,2.5] 


0.05 


-^turb 


[0.1,l.l]*5,e, 


0.1 



in the case of a constant spectral index. 



The range and binning step considered for each of the above parameters 
are given in Table [U All the other parameters of the models of the Galac- 
tic magnetic field and matter density are fixed to values proposed in Sec- 
tion m Notice that to be able to compare the dust models to the Archeops 
353 GHz data, the simulated maps are mult iplied by a mask to accou nt for 



the Archeops partial sky coverage of 30% (IMacfas-Perez et al 



2003). 



5.1. Likelihood based analysis 

To compare the Wmap data sets and the model of Galactic polarized 
emissions we used a maximum likelihood analysis where the the log-likelihood 
function is given by 

- log = 5^ {Dl^ - Ml^Nr^^,{Dl^ - M^) (18) 

i j=0 

where D^ j and M^^j correspond respectively to the data set and model for the 
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frequency of observation u. i and j index the polarization states Q and U and 
the pixel number in the maps respectively. A^~^ is the pixel-to-pixel inverse 
covariance matrices for the frequency band u. These matrices allow us to 
estimate the noise for a given pixel and the correlation between pixels. They 
are composed of 6144x6144 elements allocated in 4 blocks of 3072x3072 el- 
ements each. Each element represents the auto-correlated noise associated 



to a pixel or the correlation between 2 pixels. The e 



pixel located outer the processing mask are set to (IHinshaw et al. 



ements linked to the 



2009 



Limon et al 



20091 ). These matrices are available on the LAMBDA website]^ 



We used the processing mask of the Wmap team. This processing mask 
is built from intensity cuts on the Wmap polarized m aps at 23 GHz anc 



the thermal dust emission model of the Wmap team (IKogut et al. 



2007 



Page et al.l . 120071 ). The point-sources ha ve also been subtract e d. More de- 



tails c oncerning this mask can be found in 



Hinshaw et al. 



fl2009h 



fl2009h 



Limon et al. 



We also used a maximum likelihood analysis to compare the thermal dust 
emission model to the Archeops data. The log-likelihood function is defined 
by 



i j=0 



a2 



(19) 



where D'j^j and M^j correspond to the Archeops data and to the thermal 



^http: / /lambda. gsfc.nasa.gov/product / map / cmTent / 
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dust emission model at 353 GHz. i and j index the Stokes parameters Q 
and U and the pixel number respectively. The error bars associated to the 



Archeops data are estimated 
data (TOD) (see iPonthieu et al. 



rom 600 simulations of the 



;ime ordered 



(120051) 



Macias-Perez et al. 



mom. The 



Benoit et al 



i nstru mental noise at the TOD level was estimated following 
(120031 ). From these simulations we conclude that the noise in the Archeops 
maps can be well approximated by anisotropic white noise on the maps. We 
thus compute the variance per pixel. 



6. Results 

First of all, we consider a null turbulent magnetic field component. For 
this case, the constraints on the parameters of the polarized Galactic emission 
model, Ps, Pd and p, using the Wmap and Archeops data are presented in 
Table 121 For illustration, we also present on Figure lUthe combined marginal- 
ized likelihood at 1 and 2D using all the Wmap data. We present the 68.8% 
(black), 95.4% (dark grey) and 98% (grey) confidence level contours. The 
pitch angle, p, of the regular magnetic field component does not seem to be 
correlated neither with Ps nor Pd- Similar results are found in the case of the 
Archeops data. The best-fit values obtained for magnetic field pitch angle, p, 
are compatible for the Wmap and Archeops data, indicating that the same 
magnetic field can describe both emissions. For Pa, as it might be expected, 
constraints are only reliable at high frequency (94 and 353 GHz) and are com- 
patible within the l-cx err or bars. These results ar e also compa t ible a t la level 



with results presented in 



Boulanger et al 



(119961) 



Gold et al 



( 120091 ). Notice 



that the value obtained is for a fix degree of polarization. An independent 
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Tabic 2: Best- fit parameters for ttie models of the Galaetie polarized emissions including 
a MLS Galactic magnetic field, constrained using the Wmap and Archeops data inner 
and outer the P06 mask. The arc given by degree of freedom 



Frequency band 


zone 


p{deg) 




/3d 


Am«n 


Ka 


in P06 




-3.45 ±0.5 





10.09 


out of P06 




q p;+0.05 
'^•'-'-0.65 




^•3-0.15 


2.89 


all-sky 


-35.0^i 


_o 4P-+0.05 


0.451°:^, 


10.08 


Q 


in P06 




_o 4P-+0.05 




3.65 


out of P06 


-so.otll 


-3 fi5+°-25 

"-••"^-O.OS 


o.8i°:I 


1.58 


all-sky 




-3.45lJ]:2 




3.67 


V 


in P06 


-is.oti? 


q 4+0.15 


1 25+°-^ 


1.26 


out of P06 




— 3.95^Q 4 


1 O+0.3 


1.07 


all-sky 




q /I +0.15 
"^■^-0.85 


1 9C-+0.7 


1.26 


W 


in P06 


-60.0^?^ 


q 9+0.15 
"J-^-l.O 


-L-JO-O.IS 


1.32 


out of P06 


-35.0^^° 


q 7+0.2 
' -0.6 


9 1 c:+0.25 
^•-■-^-0.05 


1.2 


all-sky 


-60.0i?° 


q 9+0.2 
'-'•^-1.0 


1 c:+0.7 
^■'^-0.6 


1.32 


all Wmap bands 


all-sky 




-3.45l°i 




16.30 


Archeops 353 GHz 


30 % sky 


-35.01}^ 





1 O+0.4 

-L-^-o.s 


1.105 
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0,5 J.O 1.6 2.0 -U -3.6 -3.a -80 -60 -4S -20 

h f. 



Figure 1: Combined 1 and 2D marginalized likelihood, using all the Wmap channels, for 
the parameters Aturb, Ps, Pd and p assuming no turbulent magnetic field component. We 
present the 68.8% (dark), 95.4% (dark grey) and 98% (grey) confidence level contours. 
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0,5 f.d f,5 ?3 -4,5 -3,5 -3,0 -BO -60 -« -^0 




D,5 f,0 f,5 10 -t.O -3,5 -3fl -flfl -60 -« -20 0,0 0,3 04 06 06 f,0 



Figure 2: Combined 1 and 2D marginalized likelihood, using all the Wmap channels, for 
the parameters Aturb, Ps, Pd and p including a turbulent magnetic field. We present the 
68.8% (black), 95.4% (dark grey) and 98% (grey) confidence level contours. 
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determination of the degree of polarization and of the dust spectral index is 
not possible given the available data. Indeed, both parameters act as mul- 
tiplicative factors. In order to discriminate between them we would need a 
larger frequency coverage and sampling. The /3s parameter is well constrained 



at the Wmap frequencies, with a best-fit value of (5^ 



3.4;°;2, using all 



the frequency bands and it is consistent across the freq uency range. This re 



suit is cons i stent at the l a with the r esults obtained by 



Gold et al. 



( 2009 



20111 ): IPage et al.l f l2007h : 



Sun et al. 



Fauvet et al. 



mm; 



(2011); 



Jansson et al. 



(120091). 



For the Wmap data we computed the likelihood for the full sky, inside and 
outside the Wmap processing mask. In the case of Ps, the results obtained 
inner and outer the mask are consistent with those from the full sky analysis 
within the error bars. Thus, we can conclude that the spectral variations are 
smaller than the error bars in the determination of 0,^. This seems to be con- 



sistent with the results from 



Gold et al 



(12009 



201 if ) which found that the 



synchrotron spectral index is relatively constant (with respect to our error 
bars) across the sky but for the Galactic plane where they found a significant 
increased. Although the uncertainties on the dust sp ectral index are much 
l arger o ur results seem to be consistent with those in iPlanck- Collab or at ion 



(j2011bl ) that favor a larger dust spectral index at high Galactic latitudes. 



The constraints on the parameters of the polarized Galactic emission model 
including a turbulent magnetic field component are presented in Table [3] 
for the Wmap and Archeops data. As above, the best-fit parameters of 
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Tabic 3: Best- fit parameters for tlie models of the Galactic polarized emissions including 
a MLS Galactic magnetic field and a turbulent magnetic field. The are given by degree 
of freedom. 



channel 


^turb 


P 


/3. 


/3d 


Am«n 


Ka 


0.3 ±0.2 




—3. 55^*^0 5 


o.4^J:° 


13.081 


Q 




_4n+25 


-3.55^°- 


0.7l°i 


4.554 


V 


< 0.9 (95.4 % CL) 


Q + 5 


-3.6l°i 


1 oc+o.e 

^-•^^-0.3 


1.530 


w 




_c;+7 

^-37 




-■-•"-0.8 


1.463 


all 


< 0.5 (95.4 % CL) 


-30^}^ 


-3.55l°| 


1 9+0.4 
-■-•^-0.3 


20.0 


Archeops 353 GHz 


< 2.25 (95.4 % CL) 


_2Q+80 





J--«-0.9 


1.98 



the model are consistent from 30 to 353 GHz but the constraints are much 
more loosy. In particular, the relative amplitude of the turbulent component, 
Afurb, is poorly constrained although the data do not seem to favour a strong 
turbulent component. We present on Figure [2] the combined marginalized 
likelihood at 1 and 2D using all the Wmap data. We present the 68.8% 
(black), 95.4% (dark grey) and 98% (grey) confidence level contours. 
Figures [3] and H] compare the Q (top) and U (bottom) Stokes parameter 
maps obtained from the best-fit parameters (right) to the Wmap Ka and 
W band maps (left), respectively. The maps are presented in jiKju units. 
At 33 GHz the synchrotron emission dominates the signal and the data is 
well represented by the model. However, at 94 GHz thermal dust emission 
dominates the signal. This can be observed both on the Q and U maps, 
although the latter is too noisy for a clear detection of thermal dust emission. 
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Figure 3: Maps (in ^Kj^j) for the polarization Q (top) and \J{bottom) Stokes parameters 
at 33 GHz for the Wmap 5 years data (left), and or the best-fit parameters model of the 
polarized foreground emissions {right) . 
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Figure 4: Maps (in jiKnj) for the polarization Q {top) and U {bottom) Stokes parameters 
at 94 GHz for the Wmap 5 years data (left), and for the best fit model of the polarized 
foreground emissions (right) . 



22 



7. Conclusions 



We proposed in tliis pape r an extended study of tlie polarized Galactic 



diffuse emissions presented in iFauvet et al.l ( 1201 ll ) . We have constructed co 



herent models of these two foreground emissions based on a 3D representation 
of the Galactic magnetic field and of the distributions of relativistic electrons 
and dust grains in the Galaxy. For the Galactic magnetic field we considered 
a modified logarithmic spiral model for the large-scale regular component, 
plus a turbulent one. The relativistic electrons and dust grains distributions 
have been modeled with exponentials peaking at the Galactic center. 

We performed a likelihood analysis to compare the available Wmap and 
Archeops data to a set of models obtained by varying the pitch angle of the 
regular magnetic field, the relative amplitude of the turbulent magnetic field 
as well as the extrapolation spectral indices for the synchrotron and thermal 
dust emissions. From this analysis, we observe that the best-fit parame- 
ters are compatible across the frequency range explored, indicating that the 
polarized sky emission at the different frequencies are of the same nature. 
Using the full data set, we have been able to set constraints on the pitch 
ang le, p = -30t\r de grees. The best-fit value for the spectral index of the 
synchrotron emission, /3s = — 3.45^*^05, is low er but compat ible with other val- 



ues fo und in the literature (see for example 



Kogut et al 



feooTI); 



Gold et al 



(120091 )). This low spectral index indicates that the synchrotron emission ob- 
served at microwave frequency is produced by relativistic electrons with an 
steep energy spectrum (spectral index p = A). An upper limit on the relative 
amplitude of the turbulent component is obtained although it seems that 
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this turbulent part is not required to reproduce the microwave data at large 
angular scales. However, we only accounted here for a statistically isotropic 
turbulent component and we did not for the ordered turbulent component 
which in our case can not be distinguished from the regular component. 



From above we can conclude that a simple model of the polarized Galac- 
tic diffuse foreground emission including synchrotron and thermal dust can 
account for the observed sky emission in the frequency range from 30 to 
353 GHz at large angular scales. This is of great interest for the analysis 
of Planck satellite mission data, as Planck will measure the polarized 
CMB anisotropies in the frequency range from 70 to 217 GHz an d we expect 



the foreground emi s sions to be dominant at large angular scales (iGold et al. 



2011 



Fauvet et al. 



2OI2I). 
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